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Resumen 


Se analizan dos modelos mecanicistas para simular el funcionamiento hidráulico de 
un sistema de drenaje agrícola subterráneo; el primer modelo resuelve la ecuación 
de Richards sobre un dominio bidimensional, usando una condición de frontera 
tipo radiación no lineal en el perímetro del dren; el segundo modelo aplica la 
ecuación de Boussinesq con coeficiente de almacenamiento variable en un dominio 
unidimensional, imponiendo en los drenes condiciones de frontera de radiación 
fractal. Se evalúa la capacidad de descripción de ambos modelos de simulación, 
considerando el caso de un acuífero libre somero con recarga vertical nula, 
obteniéndose que la aproximación basada en la ecuación de Richards proporciona 
mejores descripciones de los procesos de transferencia de masa y energía que 
ocurren en los sistemas de drenaje agrícola, y que la ecuación de Boussinesq no 
puede describir simultáneamente las transferencias de masa y energía descritas 
con la ecuación de Richards. Se muestra, por una parte, que cuando el modelo 
de Boussinesq se usa en la resolución de problemas inversos para estimar la 
conductividad hidráulica a saturación del suelo (K) se obtienen valores para 
este parámetro mayores a los calculados con el modelo basado en la ecuación de 
Richards. Por otra parte, se denota que cuando K_ es estimada a partir del modelo de 
Richards y el valor así obtenido es usado en el modelo de Boussinesq para estimar 
la separación entre drenes, se obtienen separaciones más cortas de lo necesario. 


Palabras clave: flujo de drenaje, radiación no lineal, radiación fractal, coeficiente de 
almacenamiento variable, abatimiento de la superficie libre. 


Introducción 


Una de las funciones básicas de los sistemas 
de drenaje agrícola subterráneos es evitar el 
establecimiento de un régimen de humedad 
adverso para el desarrollo de los cultivos, 
mediante el abatimiento de mantos freáticos 
someros y el desalojo oportuno de los excesos 
de agua provenientes de sobre-riegos, preci- 
pitaciones, pérdidas por infiltración en canales 


y aportaciones de corrientes subterráneas. 
Para realizar una evacuación eficiente del 
agua es indispensable conocer cómo se mueve 
en el suelo, ya que esto permite determinar la 
evolución del caudal de agua que puede ser 
removido del medio poroso mediante drenaje y 
la variación correspondiente del manto freático. 

La determinación de las variables de 
un sistema de drenaje agrícola requiere del 
análisis de las transferencias de masa y energía 
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que ocurren en el suelo; el estudio de estas 
transferencias, aunque es altamente complejo, 
dado que se trata de analizar procesos 
básicamente no lineales que ocurren en un 
medio cuyas propiedades varían en el tiempo 
y en el espacio, puede realizarse considerando 
alguna de las tres escalas de análisis siguientes: 


Escala microscópica. En esta escala de análisis, 
que corresponde a cada poro del suelo, la 
velocidad media del agua o gasto microscópico 
unitario se estima con la ecuación de Poiseuille, 
la cual proviene de las ecuaciones de Navier- 
Stokes, y la presión en cada poro se estima con 
la ecuación de Laplace. Esta escala de análisis 
es recomendable para una comprensión fina de 
los mecanismos fundamentales de los procesos 
de transferencia. 


Escala macroscópica. La complejidad que plantea 
la definición precisa de la forma geométrica 
del espacio poral hace que la descripción 
microscópica no pueda ser puesta en práctica 
sin un cambio de escala, cuya etapa esencial 
consiste en introducir el concepto de volumen 
elemental representativo, el cual permite 
establecer una equivalencia entre el medio 
real disperso y un medio ficticio continuo. En 
esta escala de análisis, que corresponde a un 
conjunto de poros de amplia gama de tamaños, 
la velocidad media en los poros llenos con agua 
O gasto unitario macroscópico se calcula con la 
ley de Darcy-Buckingham (1907) y la presión 
del agua asociada al conjunto se calcula con 
la propia ley de Laplace aplicada al tamaño 
de poro mayor; la ecuación de transferencia 
correspondiente es conocida como ecuación de 
Richards (1931). 


Escala megascópica. En esta escala de análisis, 
que corresponde a un conjunto de suelos, la 
velocidad media o gasto unitario megascópico 
se estima con la ley de Darcy promediada, 
considerando la hipótesis de Dupuit-For- 
cheimer, concerniente a una distribución 
hidrostática de presiones, y la presión del agua 


es la proporcionada por los piezómetros; la 


ecuación de transferencia correspondiente es 
conocida como ecuación de Boussinesq del 
drenaje agrícola. 

En la práctica, el análisis del drenaje agrícola 
es realizado con la ecuación de Richards (p.e., 
Zaradny y Feddes, 1979; Fipps y Skaags, 1986; 
Saucedo et al., 2002) o con la ecuación de 
Boussinesq del drenaje agrícola (p.e., Dumm, 
1954; Pandey, et al., 1992; Gupta et al., 1994; 
Samani et al., 2007). La primera ecuación 
diferencial permite realizar descripciones de 
los procesos de transferencia que ocurren en 
las zonas saturada y no saturada del suelo; 
sin embargo, su aplicación a la escala de un 
distrito de riego e incluso de la parcela se 
ve limitada tanto por la dificultad y el costo 
del 
representar las características hidrodinámicas 


trabajo experimental requerido para 
del suelo (curva de retención de humedad y 
curva de conductividad hidráulica) como por 
el esfuerzo de cómputo necesario para modelar 
el movimiento tridimensional del agua en el 
suelo. Estas limitantes han originado que el 
estudio de los procesos de transferencia de 
agua en los sistemas de drenaje agrícola sea 
principalmente realizado con la ecuación de 
Boussinesq, aproximación que si bien considera 
de manera simplificada las transferencias que 
ocurren en la zona no saturada del suelo, es 
una herramienta útil para hacer descripciones 
aproximadas del flujo del agua en el espesor 
saturado del medio poroso. 

Recientemente se han presentado dos 
modelos mecanicistas para el drenaje agrícola, 
que mejoran las hipótesis tradicionales de 
los modelos reportados en la literatura. 
Por una parte, Zavala et al. (2005) han de- 
sarrollado un modelo para el drenaje agrícola 
subterráneo basado en la ecuación de Richards 
ecuación 


bidimensional, esta 


diferencial en la frontera del dren a una 


sujetando 


condición de radiación no lineal, con la cual 
las transferencias de masa y energía en un 
sistema de drenaje son mejor representadas. 
Por otra parte, Zavala et al. (2004) y Fuentes 
et al. (2009) han estudiado el drenaje agrícola 
con la ecuación de Boussinesq y deducido, 
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respectivamente, el tipo de condición de 
frontera que debe usarse en los drenes agrícolas 
con esta ecuación, y la relación entre la curva 
de retención de humedad y el coeficiente de 
almacenamiento en acuíferos libres someros. 

El objetivo de este trabajo es aplicar los dos 
modelos ya descritos para evaluar y comparar 
las transferencias de masa y energía obtenidas 
con las ecuaciones de Richards y Boussinesq, 
así como identificar el efecto de las hipótesis 
consideradas en la derivación de la ecuación 
de Boussinesq sobre las descripciones que 
proporciona esta relación. 


Materiales y métodos 
Escala macroscópica 


El movimiento del agua en el suelo, y en 
consecuencia, en un sistema de drenaje 
subterráneo, es un fenómeno tridimensional, 
motivo por el cual su descripción debe 
realizarse usando la ecuación de Richards en 
tres dimensiones; sin embargo, dado el esfuerzo 
de cómputo que implica resolver esta forma de 
la ecuación de Richards, es conveniente aceptar 
la hipótesis de que el fenómeno es básicamente 
bidimensional, es decir, que se efectúa en 
planos perpendiculares a la dirección del 
dren. Si además se asume que la extracción de 


agua por las plantas es despreciable, la forma 
bidimensional de la ecuación de Richards 
puede escribirse: 


1) o) 
C(w) E y 0X 


O 


(1) 


donde y es el potencial de presión del agua 
en el suelo [L]; C(y)=d0(y)/dwy, la capacidad 
específica [L”]; O(y), el contenido volumétrico 
[L9L*]; K(yw), la conductividad 
hidráulica del suelo [LT]; x, la coordenada 
horizontal 


de agua 


[El; z, la coordenada vertical 
orientada positivamente hacia abajo [L], y t es 
el tiempo [T]. 

Si se considera un sistema de drenaje con 
tubos paralelos equidistantes instalados a igual 
profundidad, se puede definir un dominio para 
la ecuación (1), como el mostrado en la figura 1. 

La descripción del drenaje agrícola sub- 
terráneo con la ecuación (1) requiere la 
inicial de las 
presiones en el medio poroso como de las 


definición tanto del estado 


condiciones en la frontera del sistema. El 
estado inicial de la presión del agua en el suelo 
se especifica como una función conocida del 
espacio: 


y=v;(x,2) (2) 


Dren (Q) 


Superficie del suelo 


Estrato impermeable 


== == O 


Figura 1. Dominio de solución para la ecuación de Richards bidimensional. 
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Por la simetría del flujo se sabe que el flujo 
de Darcy en dirección perpendicular a los 
segmentos AF, DE y BC es nulo (condición tipo 
Neumann), y una situación similar ocurre en 
el segmento de frontera CD debido al estrato 
impermeable: 


-K(w) 3-0 x=0  2€AF y zEDE 
t>0 (3) 
K(w) 20 x=1/2  zEBC 

t>0 (4) 
kw 2020 xECD  2=P 

t>0 (5) 


donde P, es la profundidad del estrato 
impermeable, medida a partir de la superficie 
del suelo; y L es la separación entre drenes. 

En la superficie del suelo (segmento AB), 
si se conoce la evolución de la evaporación (1), 
puede imponerse una condición de frontera 
tipo Neumamn: 


t>0 (6) 


De acuerdo con Zavala et al. (2005), la 
transferencia de agua del suelo al dren debe 
ser descrita con la siguiente condición de 
frontera de radiación no lineal: 


donde 2 es la derivada en la dirección normal 


a la osea del dren; q,, un valor particular del 
flujo del agua en la interfaz suelo-dren [LT”]; a. 
y [, los parámetros de forma adimensionales; 
P, la profundidad del dren [L], y h, es la presión 
en el interior del dren [L] igual a la atmosférica 
(h, = 0) en el segmento del perímetro interno del 
dren en contacto con aire, e igual al tirante de 
agua que se tenga en cada punto del perímetro 
interno del dren en contacto con agua. 

La aplicación de la relación (7) requiere del 
conocimiento de la evolución del tirante medio 
en el dren;, si se desprecian las variaciones en el 
sentido longitudinal, el tirante en el dren puede 
suponerse uniforme en el espacio, pero variable 
en el tiempo. En esta situación, la evolución en 
el tiempo del tirante puede calcularse a partir 
de una ecuación que relacione la velocidad 
del flujo con la pérdida de energía en la 
dirección del movimiento. Se considera, por 
su generalidad, la ley de resistencia fractal 
propuesta por Fuentes et al. (2004): 


e 3d-1 5d 


donde V es la velocidad media del agua en el 
dren [LT"]; d, un parámetro adimensional que 
varía entre 1/2 < d < 1, en función del tipo de 
flujo (turbulento o laminar); k, un coeficiente 
adimensional; g, la aceleración gravitacional 
[LT?]; v, la viscosidad cinemática del agua 
[L?T“]; R,, el radio hidráulico [L], y J es la 
pendiente de fricción, que en el caso de flujo 
uniforme es paralela a la pendiente topográ- 
fica del dren J, [LL]. 

La combinación de la relación (8) y la 
ecuación de continuidad para flujo permanente, 
que indica que el gasto Q es el producto del 
área hidráulica por la velocidad media del 
agua (Q = VA), permite obtener la relación 
entre el tirante medio en el dren y el gasto que 
conduce Q = kg AR39J0/v21, Sin embargo, 
la aplicación de esta relación presenta una 
limitante: se tienen dos incógnitas (gasto y 
tirante) y sólo una ecuación. Este problema se 
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resuelve planteando una segunda ecuación que 
se Obtiene al integrar la condición de radiación 
no lineal (7): 


a h B 
Q=2/ s[1-*) aa da (9) 
J Pp) UP 


donde Qes el perímetro de lasemicircunferencia 
del dren (ver figura 1) y / es la longitud de éste. 

Para modelar el drenaje agrícola con el 
sistema de ecuaciones (1-9), es indispensable 
disponer de representaciones analíticas de las 
características hidrodinámicas del suelo (y) y 
K(6). En aplicaciones de campo y de laboratorio, 
Fuentes et al. (1992) recomiendan usar el 
modelo para la característica de humedad de 
van Genuchten (1980), sujeto a la restricción de 
Burdine (1953): 


con n>2 (10) 
donde 6, es el contenido volumétrico de agua 
a saturación; 8, es el contenido volumétrico 
residual, y y,, m y n son parámetros empíricos. 

Y para la conductividad hidráulica sugieren 
el uso del modelo de Brooks y Corey (1964): 


(11) 


donde K, es la conductividad hidráulica a 
saturación [LT"*] y y es un parámetro de forma 
adimensional positivo. 


Escala megascópica 


Descripciones generales de las transferencias 
de masa y energía en los sistemas de drenaje 
agrícola pueden ser obtenidas con la ecuación 
de Boussinesq. Conforme a la hipótesis de 
que las variaciones en el sentido del dren son 


despreciables y la recarga es nula, la dinámica 
del agua en el espesor saturado del suelo en un 
acuífero libre somero puede ser descrita con la 
siguiente forma de la ecuación de Boussinesq 
del drenaje agrícola: 


(12) 


donde H y H, son, respectivamente, la elevación 
de la superficie libre y la elevación del estrato 
impermeable, medidas a partir de un nivel de 
referencia [L]; cuando el estrato impermeable es 
aproximadamente horizontal, puede suponerse 
como el nivel de referencia y tomarse H, =0; 
u(H) es el coeficiente de almacenamiento [L*L?], 
que en un acuífero libre somero es función de 
la carga hidráulica (Hilberts et al., 2005; Fuentes 
et al., 2009). 

La consideración del modelo para la curva 
de retención de humedad de van Genuchten 
sujeto a la restricción de Burdine y la hipótesis 
de una distribución hidrostática de presiones 
permite obtener la siguiente representación 
analítica para el coeficiente de almacenamiento 
(Fuentes et al., 2009): 


«060.0 01)" 


con m=1-2/n (13) 


donde H, es la elevación de la superficie del 
suelo. 

Para resolver la ecuación (12) sobre el 
dominio de solución mostrado en la figura 2, es 
indispensable definir las condiciones iniciales 
y de frontera. La especificación de estas 
condiciones límites es más conveniente si la 
posición de la superficie libre contada a partir 
del estrato impermeable es expresada como: 

H(x,t) =D, +h(x,t) (14) 
donde h(x,t) es la carga hidráulica contada a 
partir de la posición de los drenes y D, es la 
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profundidad del estrato impermeable medida 
a partir de los drenes (ver figura 2). 

En términos generales, la condición inicial 
de la presión debe ser especificada como la 
elevación de la superficie libre a lo largo de la 
coordenada horizontal x: 


h(x,0)=h(x) 


Zavala et al. (2004) han mostrado que la 


(15) 


relación que debe sujetar a la ecuación de 
Boussinesq en la frontera de los drenes es la 
siguiente condición de frontera tipo radiación 


fractal: 
2s-1 
AA (7 (=)-o (16) 


gx 
donde y es un coeficiente de conductancia 
adimensional; s, el cociente de la dimensión 
fractal del suelo D, y la dimensión del espacio 
de Euclides (s = Dj! 3); Ks, la conductividad 
hidráulica de la interfaz suelo-dren. El signo 
positivo en la ecuación (16) es tomado para 


el dren ubicado en la coordenada x = 0 y el 
negativo para x = L. De acuerdo con Zavala et 
al. (2007), es conveniente expresar la ecuación 
(16) en la forma: 


(16.1) 


donde qs = y Ks P/L. 

La dimensión cociente s está definida 
implícitamente en función de la porosidad 
volumétrica (p) del objeto fractal (suelo y pared 
del dren) como: 


(1-0) +42 =1 (17) 
Aplicación 
La comparación de las transferencias de masa 
y energía proporcionadas por los sistemas 


de ecuaciones (1-11) y (12-17) se realiza 
considerando la información de experimental 


Superficie 


. La] | E 
: Estrato impermeable ] 
xE0 x=L/2 as 


sa s 


L 


Figura 2. Dominio de solución para la ecuación de Boussinesq unidimensional. 
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(2004). El 
experimento realizado en laboratorio consistió 


de drenaje de Zavala et al. 
en lo siguiente: 


a) Construir un contenedor con hojas de 
acrílico transparente; instalar en éste, aigual 
profundidad y de forma paralela, dos tubos 
lisos de PVC con perforaciones circulares 
en su pared, para representar las ranuras 
de los drenes agrícolas; rellenar el módulo 
experimental con una muestra alterada de 
suelo de textura arenosa de la región de 
Tezoyuca, Morelos, previamente tamizada 
(malla de 2 mm) y colocada en capas 20 cm, 
manteniendo aproximadamente constante 
la densidad aparente. 

b) Tapar la descarga de los drenes y saturar el 
sistema de drenaje mediante una aplicación 
de lámina de agua en la superficie del suelo; 
dejar reposar el sistema hasta alcanzar la 
estabilidad de la presión del agua en el 
suelo; cubrir la superficie del mismo con 
un plástico para eliminar la evaporación, y 
destapar los drenes para medir el volumen 
de agua evacuado del suelo. 


Las características físicas del sistema de 
drenaje de laboratorio son L = 1.0 m, P = 1.20 
m, D, = 0.25 m; longitud y diámetro del dren, 
( =0.30 m y D,= 0.05 m, y la pendiente de los 
drenes J, = 0.001. 

La caracterización hidrodinámica del 
suelo, realizada de manera independiente 
al evento transitorio de drenaje, permitió 
determinar los siguientes valores para los 
modelos de van Genuchten, y de Brooks y 
Corey (Zavala et al., 2004): 0, = y = 0.539 cm?/ 
cm*, m = 0.373, n = 3.767, K_ = 0.183 m/h y 
y, = 0.418 m. Siguiendo la recomendación 
de Haverkamp et al. (2005), se puede reducir 
el número de parámetros a determinar en la 
caracterización hidrodinámica, asumiendo 
que el contenido volumétrico de agua es cero 
(0, =0 cm?/cm?), lo cual implica aceptar que 
cuando y —= —oo, necesariamente 8 = 0. Para 
estimar el coeficiente adimensional de la ley de 
resistencia fractal (ecuación (8)), Zavala et al. 


(2005) asimilan esta ley a la relación de Hazen- 
Williams, considerando d = 0.54, y determinan 
k = 9.83 (drenes de PVC lisos). 

Las condiciones experimentales permiten 
identificar, por una parte, que la condición 
inicial del evento de drenaje es la de un suelo 
totalmente saturado, con una distribución 
hidrostática de presiones y, por otra parte, que 
todos los segmentos de frontera del dominio, 
con excepción del asociado al dren, son 
impermeables. Introduciendo las porosidades 
volumétrica del suelo (4, = 0.539 cm*/cm?) y 
de la pared del dren (4, = 0.017 cm*/cm?*) en 
la ecuación (17), se obtiene el exponente de la 
condición de radiación fractal (16.1) para la 
ecuación de Boussinesq s= E + s,) = 0.636. 
Finalmente, los parámetros a, f$ y q, que 
intervienen en la condición de radiación 
no lineal (ecuación (7)) para la ecuación de 
Richards, así como el parámetro de escala q, de 
la condición de radiación fractal que sujeta a la 
ecuación de Boussinesq en los drenes (ecuación 
16.1), son determinados a partir de la evolución 
en el tiempo de la lámina drenada experimental. 
Para salvar el problema de estimar tres 
parámetros a partir de un solo evento, Zavala 
et al. (2005) asumen a = fB, y obtienen a = f = 
1.88 y q,=300 K. (este parámetro representa un 
promedio de la velocidad del agua en los poros 
grandes de la pared del dren y de la velocidad 
en los poros del suelo en contacto con el dren); 
en el caso de la ecuación de Boussinesq, Zavala 
et al. (2004) reportan q,= 0.913 K.. 

Para comparar las transferencias descritas 
con los dos modelos analizados es necesario 
resolver numéricamente los sistemas (1-11) 
y (12-17). Ambos sistemas de ecuaciones se 
solucionan realizando la discretización espacial 
con el método del elemento finito tipo Galerkin; 
la discretización temporal, con un método de 
diferencias finitas implícito, linealizando el 
sistema resultante con el método iterativo de 
Picard y resolviendo el sistema de ecuaciones 
algebraico con un método de gradiente 
conjugado precondicionado, adaptado para el 
manejo de matrices almacenadas vectorialmen- 
te sin coeficientes nulos. Las particularidades 
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de los métodos están documentadas, por 
ejemplo, en Zienkiewicz et al. (2005). 

La discretización del dominio de solución 
de la ecuación de Richards bidimensional se 
realiza aplicando el programa Argus One 4.0p, 
con el cual se generó una malla de elementos 
finitos de 10 795 nudos, distribuidos en 21 082 
elementos (figura 3), siendo los pasos de 
espacio mínimos Ax,  = Az, =0.2 cm y los 
máximos Ax, = Az, = 2.0 cm. El dominio 
de solución de la ecuación de Boussinesq 
unidimensional se discretizó generando una 
malla de 201 nudos y 200 elementos finitos de 
tamaño uniforme (Ax = 0.5 cm). La densidad de 
la discretización espacial de ambos dominios 
de solución permite tener errores de balance 
de masa menores del 0.1% en las simulaciones 
desarrolladas. 

Aplicando las soluciones numéricas de 
los sistemas (1-11) y (12-17), se modela el 
experimento de drenaje para determinar y 
comparar la evolución en el tiempo de la 
lámina de agua evacuada por el dren y la 
variación correspondiente de la superficie 
libre a la mitad de separación entre drenes. 


Los resultados obtenidos para la transferencia 
de masa se presentan en la figura 4 y para la 
transferencia de energía en la figura 5; en 
ambos casos, conforme el tiempo de drenaje 
aumenta las evoluciones calculadas, tienden a 
un valor límite debido a que la recarga es nula. 

En la figura 4 se observa un buen acuerdo 
entre la evolución de la lámina de agua 
drenada descrita con la ecuación de Richards 
y la evolución obtenida con la ecuación de 
Boussinesq (R? = 0.9847), resultado lógico, 
porque ambas evoluciones son el producto 
directo de la estimación de los parámetros 
de las condiciones de radiación (7) y (16.1), 
a partir de la serie experimental de la lámina 
drenada. Sin embargo, la evolución del manto 
freático no fue considerada para optimizar 
parámetros, por lo que esta variable es un 
buen elemento de comparación. Observando 
los resultados mostrados en la figura 5, 
se identifica una discrepancia importante 
entre la evolución descrita con el modelo 
basado en la ecuación de Richards respecto 
de la evolución obtenida con el modelo de 
Boussinesq; en todo el tiempo de drenaje, el 
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Figura 3. Malla de elementos finitos para la ecuación de Richards bidimensional. 
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Figura 5. Comparación de la evolución de la superficie libre a la mitad de separación entre drenes descritas 
con los modelos (1-11) y (12-17). 


abatimiento de la superficie libre calculado 
con Boussinesq es más lento que el obtenido 
con Richards. 

Tomando en cuenta que la ecuación de 
Boussinesq puede obtenerse integrando la 
ecuación de Richards en la vertical con la 
hipótesis de una distribución hidrostática 
de presiones (Bear, 1972), supuesto que no 
se cumple en las regiones cercanas al dren, 
se tiene que la descripción más precisa de la 
variables hidráulicas en un sistema de drenaje 
subterráneo corresponde a la proporcionada 


por la ecuación de Richards. Considerando 
además los resultados de las simulaciones, se 
identifica que el modelo basado en la ecuación 
de Boussinesq (ecuaciones (12) a (17)) no 
puede reproducir de manera simultánea las 
evoluciones de las transferencias de masa y 
energía descritas con el modelo basado en la 
ecuación de Richards; es decir, si reproduce 
la evolución de la masa no es factible que 
reproduzca la evolución de la energía, o si 
reproduce la evolución de la energía no puede 
reproducir la evolución de la masa. 
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La limitante de una descripción simultánea 
precisa de las transferencias de masa y 
energía con la ecuación de Boussinesq debe 
tenerse presente cuando se usa para estimar 
propiedades hidráulicas del suelo o cuando 
se aplica al diseño de sistemas de drenaje. 
Por una parte, es tradicional considerar esta 
ecuación ¡para estimar la conductividad 
hidráulica a saturación del suelo a partir de 
mediciones de abatimiento o recuperación 
del manto freático; si este es el caso, el valor 
determinado será mayor al valor real de 
conductividad hidráulica del medio poroso, 
pues de acuerdo con los resultados de este 
trabajo, la ecuación de Boussinesq describe 
abatimientos menores de los que se presentan 
en un sistema de drenaje. Por otra parte, si 
la conductividad hidráulica a saturación del 
suelo ha sido estimada a partir de pruebas de 
campo o laboratorio que consideran relaciones 
más precisas que la ecuación de Boussinesq, 
y se usa este valor de K,, en combinación con 
la ecuación de Boussinesq para calcular la 
separación entre drenes, se pueden obtener 
separaciones menores de lo que realmente 
se necesita para satisfacer el abatimiento del 
manto freático. 

Para ilustrar estas dos situaciones se 
retienen las descripciones de la prueba de 
drenaje de laboratorio obtenidas con el 
modelo de Richards y se realizan simulaciones 
numéricas con el modelo basado en la ecuación 
de Boussinesq. Se determina, primero, el valor 
de conductividad hidráulica a saturación, 
que permite aproximar con el modelo basado 
en la ecuación de Boussinesq el abatimiento 
del manto freático a la mitad de separación 
entre drenes descrito con la ecuación de 
Richards. Posteriormente se retoma el valor 
de conductividad hidráulica a saturación 
determinado en laboratorio por Zavala et al. 
(2004), y se calcula la separación entre drenes 
que requiere la ecuación de Boussinesq, para 
minimizar las diferencias entre el abatimiento 
del manto freático que describe a la mitad 
de separación entre drenes y el abatimiento 
obtenido con la ecuación de Richards. 


Los resultados obtenidos para el primer caso 
se muestran en las figuras 6a y 6b. En la primer 
figura se presenta la mejor aproximación 
de la ecuación de Boussinesq a la ecuación 
de Richards en cuanto a la evolución de la 
superficie libre a la mitad de separación entre 
drenes (R? = 0.9761), obtenida con un valor de 
conductividad hidráulica a saturación para el 
modelo de Boussinesq de K_= 0.500 m/h; este 
valor es más del doble del valor determinado 
en laboratorio (K, = 0.183 m/h), con lo cual se 
muestra que las hipótesis consideradas en la 
derivación de la ecuación de Boussinesq afectan 
de forma sensible su capacidad de estimación 
de K.. La sobrestimación del valor de K_ se 
traduce necesariamente en una sobrestimación 
de la evolución de la lámina drenada, como se 
muestra en la figura 6b. 

Para el caso del cálculo de separación entre 
drenes, se obtuvo que la mejor aproximación 
entre el abatimiento descrito con el modelo 
basado en la ecuación de Boussinesq (ecuacio- 
nes (12) a (17)) y el descrito en la ecuación de 
Richards (ecuaciones (1) a (11)), se tiene cuando 
la separación entre drenes para el primer 
modelo es 0.40 m (R? = 0.9899). La diferencia 
entre la separación real del sistema de drenaje 
experimental considerado en el modelo de 
Richards (L = 1.0 m) y la separación teórica 
requerida con el modelo de Boussinesq para 
reproducir el abatimiento fijado es 153%. 


Conclusiones 


Haciendo uso de información experimental de 
drenaje obtenida en laboratorio para el caso 
de un acuífero libre con recarga vertical nula, 
se ha evaluado la capacidad de descripción 
de dos modelos mecanicistas para simular las 
transferencias de masa y energía en sistemas de 
drenaje agrícola subterráneo, uno basado en la 
ecuación de Richards bidimensional y otro en 
la ecuación de Boussinesq del drenaje agrícola. 
Se ha determinado que con la ecuación de 
Boussinesq no se pueden reproducir de manera 
simultánea las transferencias de masa y energía 
que proporciona la ecuación de Richards; por 
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Figura 6. Comparación de los resultados obtenidos con el modelo basado en la ecuación de Richards (ecuaciones (1) a (11) 
con un K_= 0.183 m/h y el modelo basado en la ecuación de Boussinesq (ecuaciones (12) a (17)) con K_ = 0.500 m/h. 


una parte, si el modelo de Boussinesq se usa y el abatimiento del manto freático, para 
para reproducir la evolución de la masa descrita distribuir proporcionalmente en éstos el 
con el modelo de Richards, necesariamente la efecto de las hipótesis consideradas en la 
ecuación de Boussinesq describe una evolución derivación de la ecuación de Boussinesq y 
de la energía más lenta que la proporcionada obtener adecuadas descripciones simultáneas 
por la ecuación de Richards; por otra parte, si de masa y energía. Además, el diseño de un 
se reproduce con el modelo de Boussinesq la sistema de drenaje debe ser consistente; si se 
variación de la energía descrita con el modelo estima K_ con la ecuación de Richards, debe 
de Richards, la ecuación de Boussinesq estimarse la separación entre drenes con esta 
sobrepredice la evolución de la masa descrita misma ecuación diferencial, y si se estima K, 
por la ecuación de Richards. con la ecuación Boussinesq, el cálculo de la 
Los límites de la ecuación de Boussinesq separación entre drenes debe ser realizado 
se traducen en una sobrestimación de la con- con la ecuación de Boussinesq. 
ductividad hidráulica a saturación cuando Los dos modelos presentados tienen un 
esta ecuación diferencial es considerada en enfoque determinista y se aplicaron para 
la caracterización hidrodinámica de suelos, reproducir datos obtenidos bajo condiciones 
o bien en separación entre drenes menores bien establecidas en laboratorio. Sin embargo, Z 
de lo realmente necesario si la conductividad ambos modelos también pueden ser aplicados 3 
hidráulica a saturación es una dato estimado a al diseño y análisis de sistemas de drenaje 3 
partir de la ecuación de Richards. parcelarios instalados en el campo, donde una E 
Habida cuenta de las ventajas y el amplio fuerte variabilidad espacial de las propieda- E 
uso de la ecuación de Boussinesq en el des hidrodinámicas del suelo puede estar E 
estudio del drenaje agrícola, y conocidos sus presente; en este escenario, la parte esencial E 
límites de descripción, se recomienda que la antes de aplicar los modelos es realizar una E 
estimación de la conductividad hidráulica caracterización hidrodinámica adecuada del 3 
a saturación con esta ecuación diferencial se medio, ya sea mediante métodos directos o E 
realice considerando de forma simultánea inversos, que permitan la determinación de E 
en la optimización de los parámetros K. parámetros globales o efectivos que satisfagan a 
y Wy la evolución de la lámina drenada las hipótesis de los modelos presentados. 3 
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Abstract 


ZAVALA, M., SAUCEDO, H. € FUENTES, C. Description of mass and energy transfers 
in agricultural drainage systems with the Richards and Boussinesq equations. Water 
Technology and Sciences (in Spanish). Vol. 1, No. 4, October-December, 2012, pp. 27-39. 


Two mechanistic models to simulate water flow in a subsurface drainage system are analyzed. 
The first model solves the Richards equation on a tiwo-dimensional domain using a nonlinear 
radiation condition for the perimeter of the drains. The second model uses the Boussinesq 
equation with the variable storage coefficient on a one-dimensional domain, applying fractal 
radiation boundary conditions to the drains. The descriptive capacity of both simulation 
models was evaluated, considering the case ofa free shallow aquifer with null vertical recharge. 
The results indicate that the approach using the Richards equation provides better descriptions 
of mass and energy transfer processes in agricultural drainage systems and the Boussinesq 
equation cannot simultaneously describe the mass and energy transfers described with Richards 
equation. It is shown that when the Boussinesq model is used to solve inverse problems to 
estimate hydraulic conductivity at soil saturation (K), the values obtained are greater than 
those calculated with the Richards equation model. In addition, when K_ is calculated using 
the Richards model and the value obtained is used in the Boussinesq model to calculate the 
separation between drains, separations are obtained that are shorter than necessary. 


Keywords: drainage flow, nonlinear radiation, fractal radiation, variable storage coefficient, 
water table drawdotwn. 
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